Quantum statistics in complex networks 
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In this work we discuss the symmetric construction of bosonic and fermionic networks and we 
present a case of a network showing a mixed quantum statistics. This model takes into account 
the different nature of nodes, described by a random parameter that we call energy, and includes 
rewiring of the links. The system described by the mixed statistics is an inhomogeneous system 
formed by two class of nodes. In fact there is a threshold energy e s such that nodes with lower 
energy (e < e s ) increase their connectivity while nodes with higher energy (e > e s ) decrease their 
connectivity in time. 



I. INTRODUCTION 
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Recently, pushed by the need to fit the available ex- 
perimental data on a large variety of networks, statistical 
physics is addressing its attention to complex networks 
0,^|J^] and in particular to scale-free networks character- 
ized by power-law connectivity distribution. The topo- 
logical properties of these networks are related to their 
dynamic evolution and play a key role in collective phe- 
nomena of complex systems |^|,^|6j. Consequently there 
is an urgent need of a general formalism able to make a 
distinction between networks. Different approaches have 
already been proposed for equilibrium graphs Q,^] . 

In this paper we will restrict our study to inhomoge- 
neous growing networks with different quality of nodes, 
described by quantum statistics. In fact we have re- 
cently presented a growing scale-free network with differ- 
ent qualities of the nodes and a thermal noise that is de- 
scribed by Bose statistics Q . On the other hand we have 
found that a growing Cayley-tree with different qualities 
of the nodes and a thermal noise is described by Fermi 
statistics In order to be synthetic in the follow- 

ing we will refer to these two networks as the bosonic and 
the fermionic networks respectively. Given the fact that 
the solution of the dichotomy between Bose and Fermi 
statistics is an attractive topic discussed in many differ- 
ent contexts, from supersymmetry []l2| to quantum alge- 
bras ||l3|] , in the first part of the paper we compare the 
growth dynamics of the two networks. We find that the 
bosonic and fermionic network are obtained by continu- 
ous subsequent addition of an elementary fan-shaped unit 
attached in two opposite directions. While the appear- 
ance of a classical system described by quantum statistics 
is not completely new [[ljjjl^] this is the first example of 
the occurrence of two symmetrically constructed models 
following Bose and Fermi statistics respectively. Always 
having in mind the general problem of the Bose-Fermi 
dichotomy, in the second part of this work we provide 
a more realistic example of network in which the two 
growth processes coexist. This is obtained by rewiring 
a bosonic network. This complex inhomogeneous system 
has two classes of nodes with increasing and decreasing 
connectivity and is fully described by a mixed statistic 



II. SYMMETRIC CONSTRUCTION OF BOSONIC 
AND FERMIONIC NETWORKS 

The bosonic network J^] is a scale- free network in which 
each node has an intrinsic quality e from a time indepen- 
dent distribution p(e). At each timestep a new node is 
added to the network attaching m links preferentially to 
more connected low energy nodes. The probability Hi 
that a new link is attached to a node of energy £j and 
connectivity fej is given by 



n? 



(1) 



The fermionic network flO[| is a growing Cayley tree of 
coordination number m+l in which nodes have an intrin- 
sic quality e from a time independent distribution p(e). 
Nodes are distinguished between nodes at the interface 
(with connectivity one) and nodes in the bulk (with con- 
nectivity m+l). At each time step a node at the interface 
can grow giving rise to m new nodes. The probability 
that a node i grows is given by the probability pi that 
the node is at the interface (its survivability), times e@ ei 



(2) 



The dynamic of the two networks is parameterized by 
that is a characteristic of the network growth and plays 
the role of the inverse temperature, i.e. = 1/T. For 
T = the dynamics became extremal and Uf , Uf are 
different form zero only for the lowest and the highest 
energy node of the network respectively. As the tem- 
perature increases, the dynamic involves also the other 
nodes and in the T — > oo limit lif and Uf do not depend 
anymore on the energy of the nodes. 

A generic bosonic network following (^) and a generic 
fermionic network following (Q) can be constructed by 
attaching a fixed elementary unit to a number of nodes 
growing linearly with the size of the network N. 

The fixed elementary unit playing the role of the 'uni- 
tary cell' in crystal lattices, is a fan-shaped element con- 
stituted by a vertex node connected to m other nodes. 
But the way in which this unit is attached is symmetric 
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in the two networks. In the bosonic network the vertex 
of the fan is a new node attached by m links to m of the 
N existing nodes of the network. On the contrary, in the 
fermionic network the elementary unit is reversed and the 
vertex is one of the (1 — l/m)N nodes at the interface 
while the m nodes attached to it are new nodes of the 
network. Consequently both networks are constructed 
by the addition of the same elementary unit attached in 
the two opposite directions. 

The mean-field equation for the bosonic and fermionic 
network describe respectively the evolution of the con- 
nectivity ki and the survivability of the nodes. In the 
bosonic network, since every new link is attached to node 
i with probability ([!]), and m new link are attached at 
each timestep, the mean field equation for the connectiv- 
ity ki is given by 



dkj 
dt 



-PHI*. 



(3) 



where • e kj is the normalization sum of the prob- 
ability Tlf, Eq. ([!]). Symmetrically, in a fermionic net- 
work every node grows with probability Ilf given by (|^). 
Consequently the probability pi(t) that a node i is at 
the interface decreases in time following the mean field 
equation 



dpi 
dt 



Pj 



(4) 



where the denominator sum is needed in order to normal- 
ize the probability Tlf . In both networks the resulting 
structure optimizes the system by minimizing the 'free 
energy' of each node of the network €j— T log(fcj) (bosonic 
network) or ej — T\ log /o^ (fermionic network). In the two 
networks this optimization is achieved in different ways. 
In the bosonic network the low-energy nodes are more 
likely to be awarded a new link while in the fermionic 
network high-energy nodes are more likely to be removed 
from the interface. While geometrically the two networks 
are related by the reversal of the elementary unit, the 
mean field equations (^) and (^), in the case m = 1, are 
symmetric under time reversal (t — > — t) and the change 
of sign of the energies (e; — ► — ej). Self-consistent calcu- 
lations ||[hJ show that the connectivity k(t\e,t') (the 
survivability p(t\e,t'j) of a node of energy e added to the 
network at time t' , follows a power-law in time with an 
exponent dependent on its energy, 

k(t, \e,t') = m (^j fBi ) with f B (e) = e~^^\ 

p(t,\e,t')= (jY Ft) withf F (e)=e^-^. (5) 



The time reversal of the two mean-field solution implies 
here that the connectivity of the nodes always increases 



in time in the bosonic networks while the survivability of 
the nodes always decreases in time in the fermionic net- 
work. The dynamic described by (^) depends on the two 
constant p B and pp given, respectively, by the solutions 
of the two equations 



1 = f dep(e) g^ Fr ^ = J dep(e)n B (e), 
1 - — = f de P( e ) e tn<-l F ) +1 = I dep(e)n F (e) 7 (6) 

772 J 

where ns/ni?(e) indicates the bosonic and fermionic oc- 
cupation numbers respectively. Thus the evolution of 
each node of the network is completely determined by 
a number, p B or pp, defined as the chemical potentials 
of a bosonic or fermionic system with specific volumes 
v B = 1 and vf — 1 + l/(m — 1) respectively. 

The quantum occupation numbers Ufl(e) and rif (e) 
appears spontaneously in the solution of the mean-field 
equations © and (|) and assume a clear meaning when 
we look at the static picture of the networks. In fact, in 
the bosonic network the total number of links attached 
to nodes with energy e, N B (e) is given by 



N B (e) = mtp(e)[l + n B (e)]. 



(7) 



In the l.h.s. of Eq. (0) the first and second terms rep- 
resent the number of outgoing and incoming links con- 
nected to nodes of energy e. Similarly, in the fermionic 
network, the total number of nodes with energy e found 
below the interface, Np(e) is given by the difference be- 
tween all the nodes of the network and those that are at 
the interface, i.e. 



N F (e) = mtp(e)[l - n F (e)], 



(8) 



Nevertheless n B (e) and nj?(e) acquire also a very spe- 
cific role in the single time evolution of the network. In 
fact, at time t, the probability tt^ (e) of attaching a new 
link to a generic node of energy e (bosonic network) and 
the probability it^p (e) that a generic node with energy e 
will grow in the fermionic network, are given by 

(e) = £ dt'S(e - e t , ) _> p ( e )n B (e) 

4\e) = fl dt'S(e - et ,)^%^U p(e)[l - n F (e)}. (9) 

These results explain the interconnection between the dy- 
namic of the networks and their self-similar aspect. In 
fact, for the bosonic network we have that the probabil- 
ity for a new node to be linked to a node with energy 
e converges in time to the same limit as the density of 
existing links pointing to nodes of energy e. Similarly, for 
the fermionic network we have that the probability that 
a node with energy e is chosen to grow converge to the 
same limit than the density of nodes in the bulk. 

The occurrence of the two quantum statistics in the 
description of such networks is due to the fact that the 
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networks are growing by the continuous addiction of the 
unitary cell but they try also to minimize the energy of 
the system (by the choice of the node to which attach a 
new link in the bosonic network or by the choice of the 
growing node in the fermionic network). The stochastic 
model behind the construction of the two networks al- 
ways involves the choice of a node in between a growing 
number of nodes, but while in the Cayley tree a chosen 
node is removed from the interface and cannot be cho- 
sen any more, in a scale-free networks there is no limit 
to the number of links a node can acquire. Consequently 
the Cayley tree is described by a Fermi distribution while 
the scale- free network is described by a Bose distribution. 

The framework of quantum statistics clarify the rela- 
tion between self-organized critical processes and scale- 
free models. In fact, the fermionic network evolution 
in the T — > limit reduces to the Invasion Percolation 
dynamics on a Cayley tree |l6 17 lj], a well known self- 
organized process JTgj ] while the bosonic network in the 
T — ► oo limit reduces to the BA model po[ for growing 
scale-free networks. 



III. MIXED STATISTICS IN SCALE-FREE 
NETWORK WITH REWIRING 

Our purpose here is to expand on the previous results 
and to discuss systems which are governed by additional 
processes on top of the simple growth discussed before. 
For example, in real networks, in addition of the appear- 
ance of new nodes, one can observe new links as well, 
or rewiring of existing links. In fact rewiring of the link 
in a scale-free network has been used to model increas- 
ing disorder in more realistic networks 
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22 1 . We show 



that the presence of such additional processes can cre- 
ate a coexistence of Fermi and Bose statistics within the 
same system. This implies that most real systems, for 
which such additional processes are present, exist in a 
mixed state, whose statistics can be described only by 
simultaneously involving both Bose and Fermi statistics. 
It is not our purpose to model any particular system at 
this point. Thus next we discuss a simple system that 
displays this mixed behavior. 

A simple example of mixed statistics is given by in- 
troducing rewiring into a bosonic network. This net- 
work is constructed iteratively in the following way: at 
each timestep a new node and m links are added to the 
network. The new node has an energy e chosen from a 
distribution p(e) and the m links connect the new node 
preferentially to well connected, low energy nodes of the 
system. As in the bosonic network without rewiring we 
assume that a new links is attached with probability 



lit oc e- 0£i k(t\ei,ti) 



(10) 



to node i arrived in the network at time ti, with energy 
&j and connectivity k(t\ei,U) at time t. Furthermore we 



assume also that at each timestep m! edges detach from 
existing nodes and are rewired to the new node. Conse- 
quently every new node will have m + m' links attached 
to it. We assume that edges connected to high energy 
nodes are more unstable, so that the probability that an 
edge connected to a node of energy a detaches from it is 
proportional to e^ £i . Consequently, the probability that 
a node i will loose a link because of the rewiring is given 
by 



117 oc e 0ei k(t\ei,ti) 



(11) 



where U is the time node i is added in the network , is 
its energy and fc(t|e,-,f,) its connectivity at time t. The 
continuous equation describing the time evolution of the 
connectivities of the nodes is given by 

dk(t\et,ti) e-^k(t\e tl U) , eP^kftfaU) 

= m— t, — n~, r — m 



dt 



(12) 

with the initial condition 



fco = k(t\ei, t) = m + m . 



(13) 



To solve (|12|) we assume that in the thermodynamic limit 
the normalization sums Zb and Zp, given by 
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z F = J2^ £i Kt\^,tj) 



(14) 



selfaverage and converge to their mean value < Zb > e 
and grow linearly in time, with the asymptotic behavior 
given by the constants [Lb and \xp, < Zp > e , 



Zb — *< Zb > c 
Zp — >< Zp > e 



m'te ptiF 



(15) 



Using (|l5|), the dynamic equation (fL2h reduces to 



dk(t\ei,ti 



I = ( e -/3(£,-MB) _ e /3(e i -/ 1F )^MB!ilM. (16) 



dt 

Consequently we have found that the time evolution of 
the connectivity fc(t|e,,£j) follows a power-law 



with 



f \ /mi«(<0 

k(t\e,t') = k [- 



(17) 



(18) 



The characteristic difference of this network from the 
bosonic scale-free network is that the connectivity of the 
nodes, due to the rewiring process, can either increase or 
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decrease in time. In fact / m j x (e) (denned in Eq. fllS|)) 
change sign at a threshold energy value 

( s = g ■ ( 19 > 

Consequently the nodes with energy e < e s increase their 
connectivity in time while nodes with energy higher then 
the threshold e s , i.e. e > e s , have a decreasing connec- 
tivity. 

After substituting k(t\ei, U) form Eq. ( p"7| ) with f m i X (e) 
given by ( |l8| ) into (|l2|) and the sum over j with an inte- 
gral , we get the self-consistent equations for the chemical 
potential [is and [ip 

e — f3(e-(i B ) 

dep(e) 



m + m' 



dep{e) 



1 — e -/3(e-A'B) 4- g^(e-Mp) 
1 _ e -j8(e-A«s) + e /3(«-MF) ' 



(20) 
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FIG. 1. Dynamical evolution of the connectivity of nodes 
with different energies. The connectivity of the nodes always 
follows a power-law, increasing or decreasing in time depend- 
ing on the energy e and the threshold value e s . 

On the same time, the distribution of edges attached 
to nodes with energy e converges to the mixed statistics 

n mix (e)p(e) = (ro + m')p(e) i - J } - 

(21) 

while the number n+(e)p(e) of the edges stochastically at- 
tached to the nodes of energy e or the number n~(e)p(e) 
of the nodes detached from nodes with energy e are given 
by 

p — j8(e— /is) 

7i+(e)p(e) = mp(e)- 



n_(e)p(e) = m'p(e) 



1 — e -p{e-HB) + g/3(e-Aip) ' 
e /3(e-/ijr) 



1 



-f3{e—fi B ) + e /3(e~^i F ) 



(22) 



respectively. The distribution n m i X (e) appears as a nat- 
ural candidate of a mixed statistics going from the [If — > 
oo limit where n m i X {e) oc 1 + ris(e) to the — * oo limit 
where n rnix (e) oc n F (e). 



We have simulated a network with m = 2 and 



= 1 



and uniform energy distribution p[e) = 1 for e G (0,1), 
with chemical potentials \ib — 0.03, fiF — 0.51 and 
e s = 0.27. In Fig. |l| we show the connectivity of the 
nodes of the network with energy values above and be- 
low the threshold e s = 0.27. The figure shows that nodes 
with energy e < e s increase their connectivity in time 
while nodes with energy e > e s decrease their connec- 
tivity in time. In Fig. |2| we report the number of links 
attached to nodes of energy e, n m i X (e) for a system size 
N = 10 4 with the data averaged over 100 runs. In the 
same figure we report also the number of nodes stochas- 
tically attached to (detached from) nodes of energy e, 
n+ (e) (n-(e)). 
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FIG. 2. The number n m ix of edges attached to the nodes 
with energy e, the number rc+(e) of the edges stochastically 
attached to the nodes with energy e, the number n_(e) of 
the nodes detached from nodes with energy e are plotted as a 
function of energy. The simulations have been obtained with 
a uniform energy distribution in the interval [0, 1]. The data 
for 10 5 timesteps are averaged over 100 runs. 

The connectivity distribution P(k) is given by the sum 
of the probabilities P(k\e) that a node with energy e has 
connectivity k. Thus, if k > m + m' we have to sum 
over all the nodes with energy lower then the threshold 



e s , while if k < 



the summation will be over the 



nodes with energies higher than the threshold, 



P(k) = 6(k-ko))l I dep(e)—\-(^- 

ko Je<e s |/mix(e)| \ «0 

+ 9(k -k)± [ dep(e)—\-(f 

ko Je>e s \j m i X {e)\ \k 0y 



-7(e) 



-7W 



(23) 
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with 



with 



7 (e) = 1 + 1/fmUe) 



7(e) > 1 for e < e s , 
7(e) < 1 for e > e s . 



(24) 



(25) 



In the limit (3 — > all the nodes of the network evolve 
in the same way with 



fmix (^) 



2m 



= A. 



(26) 



Thus, if A > every node increases its connectivity in 
time while if A < all the nodes have decreasing con- 
nectivity. In the case A = the mean field equation 
describes a system in which the connectivities remain 
constant in time. 

On the contrary in the limit (3 — > 00 the difference be- 
tween nodes with different energy is strongly enhanced. 

We have to observe that as A goes from its highest 
value A = 1/2 to negative values, the energy distribu- 
tion goes from a pure Bose distribution to a mixed dis- 
tribution with an increasing Fermi character, i.e. with 
a decreasing Fermi potential fj,p. But it is impossible to 
reach the pure Fermi statistics in this way. In fact if we 
consider the limit m = 0, the number of links in the net- 
work is not increasing in time, and the new nodes only 
acquire edges from the rewiring process. In this case the 
connectivity of the nodes decreases exponentially as 



k(t\e,U) = A;o exp(-e 



-P(e-I*r) 



(t - U)) 



with the chemical potential defined by 



We observe that in this case the network doesn't grow 
anymore and the number of edges attached to nodes with 
energy e is simply given by the Boltzmann occupation 
factor. In this case the self-consistent equation and the 
mass conservation relation are not anymore equivalent, 
the first one reducing, in the thermodynamic limit to an 
identity. For this network the probability P(k) to find a 
node with connectivity k is given by 



P(k) = 
i.e. goes like P(k) 



dep(e)koe 



/3(e- MF ) 



(29) 



IV. CONCLUSIONS 

In conclusion we have shown the symmetry between 
the fermionic and the bosonic networks emphasizing the 
role of quantum statistics. 



These two particular evolving networks are related by 
the time reversal evident in the continuum equations de- 
scribing their dynamics and in the reversed unitary unit 
by which the two network are built. This time reversal 
implies that the connectivity increases in time while the 
survivability of each node decreases in time as an energy- 
dependent power-law. The time reversal of the single 
process generate two different structures with properties 
and dynamics only described by the functionals fiB and 
fiF, at every temperature T = 1/(3. Having introduced 
these two limit simple cases and having illustrated their 
symmetry we have shown that it is possible to construct 
a new class of networks described by a mixed statistics 
that can be applied to real systems where the two differ- 
ent growth processes coexist. 
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